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We present multi-epoch observations from the Submillimeter Array (SMA) for a sample of 171 bright blazars, 
43 of which were detected by Fermi during the first three months of observations. We explore the correlation 
between their gamma-ray properties and submillimeter observations of their parsec-scale jets, with a special 
emphasis on spectral index in both bands and the variability of the synchrotron component. Subclass is de- 
termined using a combination of Fermi designation and the Candidate Gamma-Ray Blazar Survey (CGRaBS), 
resulting in 35 BL Lac objects and 136 flat-spectrum radio quasars (FSRQs) in our total sample. We calculate 
submillimeter energy spectral indices using contemporaneous observations in the 1 mm and 850 micron bands 
during the months August-October 2008. The submillimeter light curves are modeled as first-order continuous 
autoregressive processes, from which we derive characteristic timescales. Our blazar sample exhibits no differ- 
ences in submillimeter variability amplitude or characteristic timescale as a function of subclass or luminosity. 
All of the the light curves are consistent with being produced by a single process that accounts for both low 
and high states, and there is additional evidence that objects may be transitioning between blazar class during 
fiaring epochs. 



1. INTRODUCTION 

The timescales on which high-amphtude flaring 
events occur in blazars indicate that much of the en- 
ergy is being produced deep within the jet on small, 
sub-parsec scales [II H]. Understanding if/how emis- 
sion differs between blazar subclasses (i.e., BL Lacs 
objects and flat-spectrum radio quasars (FSRQs)) 
may offer important insight into the similarity be- 
tween blazars and, furthermore, can provide con- 
straints on the formation and acceleration of the jets 
themselves. 

For the synchrotron component of blazar spectra, 
the low-frequency spectral break due to synchrotron 
self-absorption moves to higher frequencies as one 
measures closer to the base of the jet Q]- This of- 
ten places the peak of the spectrum in the millime- 
ter and submillimeter bands, where the emission is 
optically-thin and originates on parsec and sub-parsec 
scales [3] , allowing direct observation of the most com- 
pact regions near the central engine. The high en- 
ergy 7-ray emission originates as a Compton process, 
typically a combination of synchrotron-self-Compton 
(SSC) and external-radiation-Compton (ERC). De- 
pending on the source properties, the synchrotron 
photons or external photons are upscattered by the 
same population of electrons that emit the millimeter 
and submillimeter spectra. Therefore the submillime- 
ter and 7-ray emission are closely linked and give the 
full information about the source emission. 

A systematic study of the submillimeter properties 
of the entire sample of Fermi blazars has yet to be con- 
ducted and is one of the primary goals of our work. We 
present here preliminary analysis of the submillimeter 
properties of Fermi blazars detected by the Submil- 



limeter Array^ (SMA) at 1mm and 850/.tm, including 
an investigation of variable behavior and the deter- 
mination of submillimeter energy spectral indices. In 
addition, we consider the connection to the observed 
7-ray indices and luminosities. 



2. SMA BLAZARS 

The Submillimeter Array (4J consists of eight 6 m 
antennas located near the summit of Mauna Kea. The 
SMA is used in a variety of baseline configurations 
and typically operates in the 1mm and 850/im win- 
dows, achieving spatial resolution as fine as 0.25" at 
850/im. The sources used as phase calibrators for the 
array are compiled in a database known as the SMA 
Calibrator List^ [5] . Essentially a collection of bright 
objects (stronger than 750 mJy at 230 GHz and 1 Jy 
at 345 GHz), these sources are monitored regularly, 
both during science observations and dedicated ob- 
serving tracks. 

To select our sample, we identified objects in the 
calibrator list that were also classified as BL Lacs or 
FSRQs by the Candidate Gamma-Ray Blazar Sur- 
vey P CGRaBS]. Of the 243 total objects in the 
calibrator hst, 171 (35 BL Lacs and 136 FSRQs) 
have positive blazar class identifications, although 
there are three sources (J0238-H66, J0428-379, and 



^The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded 
by the Smithsonian Institution and the Academia Sinica. 
http: // smal .sma. hawaii.edu / callist / callist . html 
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Figure 1; The SMA light curves for 3C 454.3. The open 
circles represent the 850/im observations, and the open 
triangles represent the 1mm observations. 



J 175 1+096) which have conflicting classifications be- 
tween Fermi and CGRaBS. Some blazars found in the 
calibrator list have been studied extensively (e.g., 3C 
279 and 3C 454.3) but the SMA blazars have not been 
studied collectively. 

Forty-four of the objects in our total blazar sample 
were detected by Fermi and can be found in the cata- 
log of LAT Bright AGN Sources (LEAS) from Abdo et 
al. [7]. J0050-094 has no redshift in either the LB AS 
catalog or CGRaBS and is not included in our study. 
Of the 43 remaining sources, 14 are BL Lac objects 
and 29 are FSRQs, with 0.03 < z < 2.19. 

We examined submillimeter light curves for all of 
the SMA blazars, with observations beginning in ap- 
proximately 2003 (see Figure 1). Typically, the 1mm 
band is much more well-sampled in comparison to the 
850m band, but visual inspection reveals that the reg- 
ularity and quality of observations vary greatly from 
source to source. Many of the objects exhibit non- 
periodic variability, either in the form of persistent, 
low-amplitude fluctuations or higher amplitude flar- 
ing behavior. 

2.1. Submillimeter Properties 

Submillimeter Luminosities. Since we are pri- 
marily concerned with comparisons to Fermi observa- 
tions, we note that only 129 of the SMA blazars (23 BL 
Lacs and 106 FSRQs) were observed by the SMA in 
either band during the three months August-October 
2008. For these objects, submillimeter luminosities 
are calculated in the standard way: 

r A 7-i2 ^obs-Pobs s 

VeL^^ = AttD^ ^ ^ ^ , (1) 

where Dl is the luminosity distance, i^obs is the fre- 
quency of the observed band, and Fobs is the average 
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Figure 2: Variability index for our sample (top: 1mm, 
bottom: 850^m), with FSRQs as the hatched 
distribution and BL Lacs as the solid distribution. There 
is no signicant difference in the class distributions in 
either band; the "tail" to the left is populated by objects 
with errors larger than the intrinsic variability. 



flux (in erg cm~^ s~'^ Hz"-'^) over the three month pe- 
riod. We adopt a lambda cold dark matter cosmology 
with values of Hq — 71 km s^^ Mpc^^, Hm — 0.27, 
and A = 0.73. 

Energy Spectral Indices. We derive submillime- 
ter spectral energy indices from observations quasi- 
simultaneous with the Fermi observations. To be con- 
sistent with the use of , we define spectral energy in- 
dex as lyF^ = v~°'^ and calculate as from the average 
of the energy spectral indices over the corresponding 
three months. We only calculate as for the 16 objects 
(8 BL Lacs and 35 FSRQs) with observations at both 
1mm and 850/zm during this time frame. 

3. VARIABILITY ANALYSIS 
3.1. Variability Index 

We roughly characterize the level of variability of 
each source using the variability index from Hovatta 
et al. [H]: 

(^max — CTF^a^) — (-Finin + Ci^niin) ^^"^ 
(^'max - CrF„„,) + (-Flnin + 0-_F,„i„) 

Figure 2 shows the distribution for the SMA blazars. 
Objects with V < are typically unsuitable for more 
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detailed variability analysis for one of two reasons: 
(1) too few data points or (2) flux measurement un- 
certainties on the order of the amplitude of observed 
variability. It is important to note that, due to dis- 
crepancies between the sampling frequency in both 
bands, the variability indices for the 850/zm band may 
be artificially depressed due to the fact that there are 
not always corresponding measurements at higher fre- 
quencies during flaring epochs. 



3.2. First-Order Continuous 
Autoregression 

We follow the method of Kelly et al. [9] , who model 
quasar optical light curves as a continuous time first- 
order autoregressive process (CAR(l)) in order to ex- 
tract characteristic time scales and the amplitude of 
flux variations. Although flaring behavior is not typi- 
cally thought of as an autoregressive process, we find 
that the light curves are well-fit by the models and 
therefore adopt the method here to study blazar sub- 
millimeter light curves. 

The CAR(l) process is described by a stochastic 
differential equation 9 , 

dS{t) ^ -S{t)dt + aVdte{t) + bdt, (3) 

T 

associated with a power spectrum of the form 



Pxif) 



2^2 



1 + (27rT/)2 ■ 



(4) 



In equations [3] and |4] r is called the "relaxation 
time" of the process S{t) and is identified by the 
break in Px{f)- The power spectrum appears flat 
for timescales longer than this and falls off as 1/ for 
timescales shorter than the characteristic timescale of 
the process. 

Taking the logarithm of the blazar light curve (in 
Jy) to be S{t), we adopt r (in days) as the character- 
istic timescale of variability, after which the physical 
process "forgets" about what has happened at time 
lags of greater than r. The two other relevant pa- 
rameters, a and fj, — b/a, are the overall amplitude 
of variability and the logarithm of mean value of the 
light curve, respectively. 

In the routine, we construct an autoregressive 
model for the light curves for a minimum of 100,000 
iterations and calculate the value of r from the break 
in the power spectrum in each instance. Due to the 
limited number of observations in the 850/im band, 
we performed this autoregressive analysis only for the 
1mm light curves, which typically have more than 10 
points per light curve. 

This method yielded some surprising results. In 
Figure 3, we see that the BL Lacs and FSRQs exhibit 
virtually no difference in characteristic timescale, with 
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Figure 3: Characteristic timescale (days) versus 
submillimeter luminosity (erg s~^) in the 1mm band for 
all objects. Physically, r represents a "relaxation 
timescale" , the timescale beyond which events are no 
longer correlated. 



both classes extending across a large range in r. Be- 
cause of the uncertainty for objects with shorter char- 
acteristic timescales, it is hard to draw any definitive 
conclusions about the differences between classes. It 
is important to note that r does not necessarily rep- 
resent a flaring timescale, which is a behavior that 
typically operates on a scale of ~10-100 days and not 
on the longer timescales we see in r. 



4. CONNECTION WITH GAMMA-RAYS 

In general, we find that in the submillimeter, we 
are observing these blazars at or near the peak of the 
synchrotron component (as ~ 0), but that Fermi- 
detected sources have more negative energy spectral 
indices overall than Jermi-nondetected sources. In 
Figure 4, we see that while the majority of Fermi 
blazars are observed on the rising part of the syn- 
chrotron component (at lower energies than the peak) , 
all of the objects have very steeply falling 7-ray energy 
spectral indexes, putting the 7-ray peak at lower en- 
ergies than the observed Fermi band. Knowing that 
we are not observing the synchrotron and 7-ray com- 
ponents at analagous points in the spectrum may al- 
low us to better understand the magnetic field in the 
parsec-scale jet region and the population of external 
photons that is being upscattered to 7-rays. 

In Figure 5, the ratio between and vL^ i^^m re- 
flects the division between BL Lacs and FSRQs as well 
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Figure 4: The 7-ray index versus submillimeter index plane. The blazars faU more steeply in the 7-rays than in the 
submillimeter band, where most are, in fact, rising. This LAT-detected sample contrasts with the full SMA sample, 
where the blazars are more distributed around qs ~ 0. 



as the presence of SSC versus ERC. Here, we use sub- 
millimeter luminosity as a proxy for jet power, which 
is correlated with the integrated luminosity of the syn- 
chrotron component. Elevated 7-ray luminosity with 
respect to the synchrotron component (which is often 
seen in FSRQs) suggests the upscattering of external 
photons off the synchrotron-emitting electrons. These 
objects should occupy the upper right of the ratio/jet 
power plot, and BL Lacs, which generally exhibit com- 
ponents with roughly comparable luminosities, should 
occupy the lower left. It is clear from the figure, how- 
ever, that many FSRQs exhibit ratios similar to those 
of the BL Lacs and vis versa. 

Sikora et al. [10] report that, during its flaring 
epochs, 3C 454.3 transitions from its typical FSRQ 
state to a more BL Lac-like state, where the syn- 
chrotron component emits much more strongly com- 
pared to the 7-ray component than during its "low 
state". 3C 454.3, which is the highest submillime- 
ter luminosity FSRQ in our sample, would then shift 
down and to the right in Figure 5 when it enters a 
flaring period. For the first three months of the Fermi 
mission, 3C 454.3 was not flaring, which may explain 
its present location in Figure 5. The three objects for 
which there is a type discrepancy between CGRaBS 
and LBAS are all FSRQs (in CGRaBS) and exhibit 



low luminosity ratios and high luminosity, which sug- 
gest they may be undergoing the same changes as 3C 
454.3. A possible interpretation of the elevated lumi- 
nosity ratios observed in some BL Lacs objects is that 
there has been a dramatic increase in 7-ray luminos- 
ity due to ERC, which would not be reflected in the 
synchrotron component. 



5. CONCLUSIONS 

The motivation for observing blazars in the sub- 
millimeter is to study behavior close to the central 
engine, where the jet material is presumably still be- 
ing accelerated. The separate emission processes that 
contribute to overall SED may present differently in 
BL Lacs and FSRQs, allowing us to understand the 
similarities and differences between blazar types. We 
have investigated these differences between objects in 
terms of submillimeter behavior and, in conclusion, 
find that 

• The SMA blazars exhibit submillimeter energy 
spectral indexes that follow the spectral se- 
quence interpretation of blazars. 
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Figure 5: Ratio of 7-ray luminosity to submillimeter luminosity in the 1mm band. The location of an object in this 
plot should be directly correlated with its blazar "state" , with FSRQs occupying the upper right and BL Lacs the lower 
left. Flat-spectrum radio quasar 3C 454.3 is the object with the highest submillimeter luminosity in this plot. 



• BL Lacs and FSRQs do not exhibit significant 
differences in amplitude of submillimeter vari- 
ability or characteristic timescale, but our sam- 
ple of BL Lacs may be dominated by high- 
peaked BL Lacs (HBLs), which exhibit obser- 
vational similarities with FSRQs. 

• Blazar submillimeter light curves are consistent 
with being produced by a single process that ac- 
counts for both high and low states, with char- 
acteristic timescales 10 < Tj-est < 500 days. 

• The blazars detected by Fermi have synchrotron 
peaks at higher frequencies, regardless of sub- 
millimeter luminosity. 

• FSRQs exhibit higher ratios of 7-ray to sub- 
millimeter luminosity than BL Lacs (Figure 5), 
but all objects inhabit a region of parameter 
space suggesting transitions between states dur- 
ing flaring epochs. 

As Fermi continues to observe fainter sources, the 
sample of objects for which we can perform this type of 
analysis will increase and provide better limits on our 
results. To understand the physical relevance of these 
results, however, it is important to be able to distin- 
guish between the difference in variability between BL 



Lacs and FSRQs. One avenue for exploring this dif- 
ference is to monitor changing submillimeter energy 
spectral index and the ratio of 7-ray to submillime- 
ter luminosity as functions of time. The full mean- 
ing of the results of our autoregressive method is not 
yet clear, and will require better-sampled blazar light 
curves and the comparison between Tjcst with physical 
timescales such as the synchrotron cooling timescale. 
These analyses would allow us to place constraints 
on the processes occurring near the base of the jet in 
blazars and further understand the intimate connec- 
tion between them. 
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